CRITICAL NONLINEAR SCHRODINGER EQUATIONS WITH 
AND WITHOUT HARMONIC POTENTIAL 



REMI CARLES 



Abstract. We use a change of variables that turns the critical nonlinear 
Schrodinger equation into the critical nonlinear Schrodinger equation with 
isotropic harmonic potential, in any space dimension. This change of variables 
is isometric on L 2 , and bijective on some time intervals. Using the known re- 
sults for the critical nonlinear Schrodinger equation, this provides information 
for the properties of Bose-Einstein condensate in space dimension one and two. 
Wc discuss in particular the wave collapse phenomenon. 



1. Introduction 

Bose-Einstein condensation is usually modeled by a nonlinear Schrodinger equa- 
tion with harmonic potential (see e.g. Q), 

m ^ + ^-Aip h = ^Lu 2 x 2 ip h + \ip n \ V\ (t, x)£lxl", 

2m 2 m 

where u> > and a is the scattering length, whose sign differs according to the 

chemical element considered. For instance, it is negative for 7 Li atoms (j||, ^|), 

as well as for 85 Rb, and positive for 87 Rb, 23 Na and 1 H. The harmonic potential 

x 2 models a magnetic field whose role is to confine the particles (this is one of the 

ingredients for Bose-Einstein condensation, once the atoms have been cooled by a 

laser, see e.g. jjj), and the nonlinear term takes the (main) interactions between 

the particles into account. To simplify the mathematical analysis, we assume from 

now on that m — h = 1, and we denote ina by A S M. Notice that in || and Q, 

we considered the semi-classical limit h — > 0. 

In the above equation, the nonlinearity is cubic, regardless of the space dimension 

n > 1. Other models are also considered. In the authors propose a quintic 

nonlinearity in space dimension one, and in [ pT| , the author suggests more generally 

the study of 

!1 u 2 x 2 
id t u + -Au = ~^-u + A|u| 4/ "u, (i, ijelxl", 
«|t=0 = u 0- 

As noticed in pH , the proposed nonlinearity is the usual critical nonlinearity for 
the nonlinear Schrodinger equation with no potential (ui — 0, see e.g. When 
n = 1, this suggestion meets the model proposed in ]ic(| , and when n — 2, this is the 
usual cubic nonlinearity. The results in pl[ enlighten a rather surprising analogy 



between the study of (1.1) and that of 
(1 , )( id t v+^Av = \\v\ 4/n v, (t,x)ERx 

«|i=0 = "0- 
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Many results are known for (1.2), we recall some of them. In pQ], the author proved 
that if 

« e2:= H\m. n ) n {4> £ L 2 (R n ); \x\<f> £ L 2 (MJ 1 )} , 

then there exists T > such that v £ C(] — T, T[, £). If A > 0, then one can take 
T = oo. When A < 0, one can take T = oo when ||Mo||i 2 is sufficiently small. More 
precisely, let Q denote the ground state, which is the unique radial solution of (see 

H, @) 

lAQ + Q = -X\Q\ 4/n Q, inR", 



(1.3) 



Q > 0, in 



Weinstein proved that if ||wo||l 2 < ||Q||l 2 j then one can take T = oo. On the 
other hand, if ||u ||l 2 > IIQIIl 2 ; then the wave v may collapse in finite time. Zhang 
proved that the same holds for the solution u of (1.1). We use a change of variables 
that shows why this is so. 



Fix co > 0. Let v be a solution of (1.2), for \t\ < T, and define, for \t\ < 
arctan(cjT , )/w, 



(1.4) 



u(t, x) 



1 



(cos cot) n l' 



tanwt 



.r 



cos cot 



Then u solves (1.1). This was first noticed in [|16| for the linear case (A = 0), and in 
|l8| for the nonlinear case with critical nonlinearity. Reciprocally, if u solves (1.1), 
then v, defined by 



(1.5) v(t,x) 



1 



(1 + {tot) 2 ) 



I arctanwi 



solves ( |l.2| ). The transforms ( |l.4| ) and (L5) do not alter the initial data uq, and 
are isometric on L 2 (W l ). Therefore, it is not surprising that global existence is 
obtained whenJuolU 2 < IIQIU 2 f° r both cases, provided that it is known for one of 
them. In Sect , g, we recall the known results for Eq. (1.2), and analyze the relations 
between (1.1) and fll.2| ) in Sect. ||. Finally, we investigate the consequences of the 
results known for (1.2) as far as Bose-Einstein equations are concerned, in Sect. 0. 



2. Known results for the critical nonlinear Schrodinger equation 

Start with the initial value problem. It is now classi cal that if uq £ iJ 1 (R"), 
then there exists a solution of the initial value problem (|l.2| ) , continuous on some 
(possibly small) time interval, with values in _ff 1 (R") (see e.g. R]). The same holds 
if uq £ £ (which is the natural space to study (1.1), since £ = D (V— A + a; 2 )). 
Moreover, the mass and energy associated to the equation are conserved. Because 
we consider the nonlinear Schrodinger equation with critical power, the pscudo- 
conformal conservation law is in our case an exact conservation law. 

Proposition 2.1. For every uq £ H 1 (M. n ) (resp. uq £ SJ, there exist T*(uo), 
T*(uo) > and there exists a unique, maximal solution 



v eC(]-r4 W o),T*( Uo )[,i/ 1 (R"))nc 1 (]-T,K),T 1 *( U o)[,^ 1 (R n )) 

(resp. v eCQ -T,( Wo ),T*K)[,S)nC ,1 (] -T,( Uo ),T*( Uo )[,H- 1 (R"))) 
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of problem (l.i). The solution v is maximal in the sense that if T*(uq) < oo, 
then \\v(t)\\jji — ► oo as t f T*(uq), and if T*(uo) < oo, then \\v(t)\\ H i — > oo 
as t 4 — T*(uo). In addition, we have the following three conservation laws for 

tG]-r*(uo),r*(«o)[. 

I. Conservation of mass: ||v(^)||l 2 = || u o||l 2 - 
Conservation of energy: 

E 1 {t) := \\\V x v(t)\\h + j^^IK*)!!^?- = ^(0). 
3. Pseudo-conformal conservation law: 

Mt) ■= \\\(x + itV x )v\\l 3 + ^-11^)11^" =^2(0). 

Remark 2.2. The lower regularity uq G L 2 could be considered as well, using the 
results of Cazenave and Weissler 

In some cases, it is known that we have T*(uo) = T*(uo) = oo (see e.g. [[?)). 

Proposition 2.3. Let uq G ff x (K") (reap, uq € Ej, and Ze< v be the maximal 



solution of (l.i ). Then v is defined globally in time, that is T*(uq) = T*(uq) = oo, 
in either of the following cases. 

• If the nonlinearity is repulsive, A > 0. 

• If the nonlinearity is attractive, A < 0, and the mass of uq is sub-critical in 



the sense that ||wo||i 2 < IIQIIl 2 ? where Q is the ground state defined by (l.l) 



On the other hand, we have some sufficient conditions for which it is known that 
the solution blows up in finite time. We restrict to the case T*{u$) < oo, which 
corresponds to a wave collapse in the future. 

Pro position 2.4 ( [po| , Th. 4.2). Recall that E\ denotes the energy associated to 



(l.i.), that it is constant. Let either 

(i) E 1 < 0, 

(ii) Ei = and Im J Hox.Vuodx < 0, 

or 

(Hi) Ei > and Im J ugx.'Vuodx < —2\/E~i\\xuo\\ ^2 . 
Then there exists < T < oo such that 

lim ||V x u(*)IU 2 = oo. 

We now assume that A < 0: global existence is not guaranteed. When blow 
up occurs (which could be the case under other conditions than those stated in 



Prop. 2.4), Merle analyzed very precisely its mechanism, when the mass is critical. 
From Prop. [2^, global existence is ensured when ||uo||l 2 < IIQIU 2 - Weinstein [ pp| 
proved that if ||uo||l 2 = ||Q||l 2 j then wave collapse may occur, at least for some 
particular initial data. Merle proved that up to the invariants of ( |l.2| ), the blowing 
up solutions enlightened by Weinstein are the only ones. 

Theorem 2.5 (Q, Th. 1). Let A < 0, u G H 1 (W l ), and assume that the solution 
v of fcLW blows up in finite time T > 0. Moreover, assume that \\uq\\l2 = \\Q\\l 2 , 



where Q is defined by (l.S). Then there exist 9 G R, S > 0, Xq,Xi G R" such that 
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and for t < T, 

v(t,x) = (^-^j n/2 e^-i\ X - Xl |V2(T- t)+ «V(T- t )Q ^ (E^l _ Xo 

Remark 2.6. In addition, Merle proved that when the mass is critical, only three 
causes can prevent the global definition of v with optimal dispersion of the L 2+4 / n - 
norm of v (@, Cor. 1.2). 

• The initial data of Th. 2.5, that cause blow up at some positive time. 

• Their conjugates, that cause blow up at some negative time. 

• The solitary waves, caused by 

u {x)=8 n ' 2 e ie Q{8{x-x Q )). 



We conclude this section by recalling another consequence of Th. 2.5. 
Corollary 2.7 ([0, Cor. 1.1). Let A < 0. The solutions of the Cauchy problem 

Ud t v + ^Av = X\v\^ n v, fort>0, 

\ \v(0,x)\ 2 = ||Q||| 2 4=o, 
are exactly, for 6 G R. S > 0, x £ R", 

/ c \ n/2 , / c 

v{t,x)=(-) e i9+i T-%Q'' 



tj -\i X - Xa 
3. Transformation and relation between the two equations 



Let uo £ E. From Prop. |2.l| , the problem (fb^) has a unique solution v £ 
CQ - T, T[, E) for some positive T. For \t\ < arctan(o/T)/w, let w be defined by 
(|L4|), that is 



(coswi)™/ 2 \ w ' cosujt / 

Then u £ CQ — arctan(q>T) /oj, arctan(cjT)/cj[, E), and u solves (IT). Co nversely, 
if u £ CQ - t, t[, E) solves (|L~i"l ) with < r < 7r/2, then v defined by ( [F5|) , that is 

/ arctanwi a; \ 



(i + M) 2 ) n/4 V w ' Vi + M) 2 / ' 

is such that w G C(] — tan(a;T)/w, tan(wr)/w[, E). Transposing similarly the results 
of Prop. 2.1 yields the following corollary, 

Corollary 3.1. Let uq £ E. Then there exist t*(uq), t*(uq) > and there exists 
a unique, maximal solution 

u £ C{) -t*(u ),t*(u )[, E) n C 1 {}-n(u ),T*{u )lH- 1 {R n )) 

of problem (1A ). It is maximal in the sense that is t*(uo) < oo, then \\u(t)\\ H i — > oo 
as t | t*(uq), and if t*(uq) < oo, then \\u(t)\\ H i — > oo as t I —t*(uq). In addition, 
the following three quantities are constant for t 6] — t»(«o)i t *(uq)[- 

1. Conservation of mass: ||u(t)||^2 = 1 1 1 1 z- 2 - 

2. First part of the energy: 

(3.1) E±(u) = ~ \\ujx sin ujtu(t) — i coswiV x u(t)\\ 2 L2 H — -j-^ cos 2 CLi^||u(t)||^{™„. 
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3. Second part of the energy: 

(3.2) Eziu) = — \\uix cos ujtu(t) + isimotV x u(t)\\ 2 L2 H — - sin 2 



Remark 3.2. The existence part of this result was proved in |17| , and revis ited in 
|6). In the special case of a critical pow er n onlinearity, the transform (1.4) shows 
that no new proof is needed when Prop. 



2.1 



is known. 



Proof. The only point that we have to prove is the blow up case. Assume for 
instance that t*(uq) is finite. Up to a time translation, we can suppose that < 
t*(uq) < n /2u> . Then v, defin ed by (1.5), solves (1.2). Let v be the maximal 
solution of (1.2) given by Prop. 2.1. If it were globally defined, then v would also 



be globally defined; the transform (1.4) would then make it possible to define u up 



to some time r such that t*(uq) < t < 7t/2oj, which contradicts the maximality of 
t*(uq). Therefore, there exists T* < oo such that 

\\^xV(t)\\ L 2 — -> OO. 

We prove that u blows up at least before time arctan^T*)/^. For < t < 
arctan(cjT , *)/cj, 



> 



—iux sin(wi)i> 
, tanwi 

VxV 



tancjf 



L- 



lux sin(cji)u 



/ tanujt 



tanwi 



L- 



L- 



Now we know that for t < T* (see for instance [E0| and references therein), 

d 2 



dt 2 



\\xv(t,x)\\ L 2 = AEi, 



where E\ denotes the energy of v, which is constant. Then letting t go to 

arctan loT* 



yields the blow up part of the corollary. 



□ 



Remark 3.3. Both conservation laws (3.1) and (3.2) were derived in |6[, in the case 
of a more general nonlinearity, not necessarily critical (the second terms of Ei(u) 
and E2{u) have to be adapted according to the power considered); in general, E\(u) 
and E2(u) do depend on time, they are constant only in the case of a critical power. 
On the other hand, the sum of E\( u) a nd E2(u) is always constant, and corresponds 
to the usual energy associated to (1.1), 



E(u) = - ||V 



\u>xu(t)\\ 



L- 



\ U (t) 11^2+4/n 



Notice that the energy for v is always conserved as well (it reflects the Hamiltonian 
structure), while the pseudo-conformal conservation law is in general an evolution 
law, which is an exact conservation law only in the critical case (and the free case 
A = 0). 
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4. Transposition of the results for Bose Einstein equations 



In this last section, we investigate some consequences of the transform (1.4) and 
its companion fll.5p for Bose Einstein equations. We consider in particular the 
questions of global existence and wave collapse phenomenon. 

Transposing Prop. 2.3 yields the following corollary, 



Corollary 4.1. Let uq £ E, and u the maximal solution of (1.1 ). Then u is defined 
globally in time, that is t^(uq) — t*(uq) — oo, in either of the following cases. 

• If the nonlinearity is repulsive, A > 0. 

• If the nonlinearity is attractive, A < 0, and the mass of uq is sub-criti cal i n 
the sense that ||wo||l 2 < ||Q||l 2 j where Q is the ground state defined by (Li). 



Remark 4.2. As mentioned in the introduction, this result was proved by Zhang 
pi"] . We believe that our approach provides a good explanation for this result. 

Proof. The repulsive case is straightforward, since the three conservations stated 
in Cor. 3.1 provide a priori bounds on the E-norm of u. When the nonlinearity 
is attractive and the mass of uo is sub-critical, we use Prop. 2.3. From the results 
of Weinstein |2(|, the solution of (1.2) is defined globally when the mass of uq is 
sub-critical, HmoIIl 2 < IIQIIl 2 - Then u £ C([0, 7r/2w[, E), from the transform (1.4). 



Recall that the L 2 -norm of u is constant on [0,7t/2cj[. Considering t$ = tt/Auj as 
a new time origin and repeating this procedure shows that u £ C7([0, 37r/4w[, E). 
Using this argument indefinitely shows that u £ C([0, oo[, E), and similarly, u £ 
CQ - oo,oo[,S). □ 
Like for the case of ( |1.2| ), we can state sufficient conditions where blow up occurs. 
This is done by transposing Prop. 2.4 with (jl.4|). 



Corollary 4.3. Recall that E\ is defined by ( RjU and is constant as long as u 
belongs to E. Let either 

(i) Ei < 0, 

(ii) Ei = and Im J UQX.\/uodx < 0, 

or 

(Hi) Ei > and Im J u^x.Wuodx < — 2v / £o'||£Mo||l 2 ■ 
Then there exists < r < ir/2(v such that 

lim || V 2: M(t)||i2 = oo. 



Remark 4.4. The above criteria yield wave collapse at time r < ir/2u>. It is sensible 
to expect this phenomenon to o ccu r possibly at time r = 7r/2w, which corresponds 
to a focus fo r th e free equation (1.1) with A = 0. This geometric aspect is hidden in 
the case of (1.2), since it corresponds to infinite times. We will consider this point 
more precisely later, in the case of a critical mass (||uo||l 2 = IIQIU 2 )- As noticed 
in HJ] and Q, wave collapse for u always occurs at time r < ir/2uj when Ei = . 
Thus we could say that the compactification of time in the transformation (1.4) 
leads to new blowing up solutions. 

On the other hand, if A < and ||wo||l 2 = IIQIIl 2 , then wave collapse occurs. 
This can be analyzed very precisely thanks to the results of Merle. As a corollary 
of Th. 2.5, with the change of variable 8' = <5coswr, we have the following, 
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Corollary 4.5. Let A < 0, uq G H 1 (R n ), and assume that the solution u of (1.1) 
blows up infinite time < r < tt/2lo. Moreover, assume that ||iio|| z, 2 = IIQIIl 2 - 
Then there exist 9 G K, 6 > 0, xq, x\ G M" suc/i i/iat 

( 4.i) Uo(a; ) = ( J±X'\ie + i^-^-^^r Q ( 5 ( *jz*i _ *<A 
\sm lot J \ \smuT to J 

and for < t < T, 

f. \ id+i8 2 " COB tt— i%\ — s_ — zi I 2 CQ . 3 — iujz2- tan uit „ 

U(t,X)=e n»»Tim«(r-l] 2 | cos ut I >m»(T-t) 2 x 



n/2 



X — X! COS Wf Xq 



sinw(T — i)y \ \sinw(r — t) uo 

Example 4.6. Blow up at time r = 7r/4o;, with critical mass |] ito II i. 2 = IIQIIl 2 j i s 
caused by initial data of the form 

u (x) = WV^^^Q (6 (x-xt- x )) . 



Recall that the quadratic oscillations always cause a focus for ( |l.2|) (see e.g. g). 
On the other hand, the geometry of the harmonic potential creates a focus at time 
it/2lo. Therefore we can say that in the above case, both phenomena cumulate, to 
anticipate the "usual" blow up. 

Example 4.7. Using a time translation shows that blow up at time r = 37r/4w, with 
critical mass ||uo||l2 = ||Q||l 2 ) is caused by initial data of the form 

u (x) = WV+^^^Q («S {x-xx- as )) ■ 

In that case, the oscillation e luJ ~ tends to delay the concentration (it is the outgoing 
oscillation once the focus has been crossed in [Q ) , but the geometry of the harmonic 
potential counterbalances this phenomenon and eventually causes wave collapse. 

Remark 4.8. We could of course state the analo gue of this result for blow up in 
the past, —it/2lo < t < 0. Since the transform ( |l.4| ) is "almost" 7r-periodic, and 
using the fact that the ground state Q is spherically symmetric, we can deduce in 
particular that if u(t) does not blow up for t G [0, tt/u>[, then it will never blow up 
in the future, and has not collapsed in the past. 



From Remark 2.6, three causes can prevent the global definition of solutions of 



(1.2) with optimal dispersion of the L 2+4 /"-norm of v, when the mass is critical. 
We have not analyzed the last possibility yet. In the case of solitary waves, we have 

v(t, x) = S n/2 e w Q (5(x - x )) e lSH . 



The transformation (1.4) then yields 

n/2 



/ 6 \ iV+iu tan ut ^Hy — ' / / 



COS LOt 



•< - 



From pH, Cor. 1.2, all the initial data uq with ||wo||l 2 — IIQIU 2 different from 



(O) and (JO) yield a solution v globally defined, with |K*)II rtw» = 0(i~ 2 ^ 



Back to u, thanks to the transformation ( |1.4| ), this provides a uniform estimate for 
ra+4?»j when t G [0,7r/2w[. From the conservation of the energy E(u) (the 
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usual Hamiltonian) , along with the conservation of mass, this yields an a priori 
estimate for the E-norm of u(t, .). Therefore, u does not blow up at time tt/2uj 
(otherwise, its i? x -norm would not be bounded near ir/2u), see e.g. 0, Th. 4.2.8). 
From Prop. fo[ there exists some positive a such that u is defined for t € [0, it/2uj + 
a[. One can use similar arguments for negative times, and using a time translation, 
we can deduce a similar description for t 6 [0,tt/oj[. Finally, if the solution has 
not blown up when reaching t = tt/oj, then it will never blow up, as mentioned in 
Remark 4.8. To summarize, we have the following, 



Corollary 4.9. Let A < 0, and uq € E be such that || iio I! i 2 — IIQIIl 2 - Assume in 
addition that 

| x — &i I 2 

~Q (S(x- xx - xq)) 



u (x) £ 5 n/2 e 



iO — iu) cot u)T - 



for (5,0,t,x o ,Xi) £ W + 
[0, 7t/2cj[. If in addition, 



7T 

°<2^ 



Then u(t) is defined for t € 



u (x)^S^ 2 e ie Q{S(x-x )) 



for (S,9,x ) e R+ x K x K n , then u(t) is defined for t S [0,7r/2w + a[ for some 
positive a. Finally, if moreover 



U 



(x) ^ ^«/2giS — iwcotur 



l*-*ir 



for {8,9,t,xq,x\) e 



0, 



(S (x — x\ 

xK"xIR™, thenu(t) is defined for t e] - oo , oo[. 



Adapting Cor. 2.7 yields the next corollary. 
Corollary 4.10. Let A < 0. The solutions of the Cauchy problem 




idtu + — Au 



X\u 



4/", 



2 

u(0,x)\ 2 = \\Q\\ L 20 X=Q , 
S > 0, x e R", |r| < tt/2uj, 

X \ ™/ 2 



/or t > 0, 



sin urf 



^+tanu(r- i)) 



and 



u(t, x) 



LOS 

sin cot 



x Q (luS 



n/2 



sin tot 



a; 



ip+iu; cot uj£ ^ — 



a; 



sin uit 



Remark 4.11. When letting w go to zero in Cor. 4.5 and Cor. 4.1C, we retrieve the 
initial results of Merle. Letting r go to k/2oj in the first part of Cor. [4. 9| yields the 
second part of Cor. 4.9, but the same does not hold in Cor. 4.10| . The two kinds of 
solutions have a different structure (it is so in the case with no potential). 

Remark 4.12. Continuation after blow-up time. In flifl , Merle con side rs the pos- 
sible continuations of the soluti on a fter the breaking time. With ( |l.4|) , we could 
adapt this theory to the case of (1.1). However, it seems very likely that ([l.]f) does 
not remain a good model for Bose-Einstein condensation after the wave collapse. 
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Remark 4.13. In (see also |l5| for other results), it is proved that one can fix k 
points in R" and construct a solution of (|l.2|) that blows up exactly at these points. 



This could be done for (1.1) as well. 
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